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We present the first study of parametric resonance in quantum field theory from a complete next-to- 
leading order calculation in a 1/iV-expansion of the 2PI effective action, which includes scattering 
and memory effects. We present a complete numerical solution for an 0(iV)-symmetric scalar theory 
and provide an approximate analytic description of the nonlinear dynamics in the entire amplification 
range. We find that the classical resonant amplification at early times is followed by a collective 
amplification regime with explosive particle production in a broad momentum range, which is not 
accessible in a leading-order calculation. 



PACS numbers: 11.15.Pg, 98.80.Cq, 05.70.Ln 

In quantum field theory the phenomenon of parametric 
resonance describes the resonant amplification of quan- 
tum fluctuations, which can be interpreted as particle 
production. It provides an important building block for 
our understanding of the (prc)heating of the early uni- 
verse after a period of inflation Q . It has been frequently 
discussed for relativistic heavy-ion collisions in the forma- 
tion of disoriented chiral condensates 0] , or the decay of 
Polyakov-loop condensates , or parity-odd bubbles . 

Despite being a basic phenomenon that can occur in a 
large variety of quantum field theories, parametric reso- 
nance is a rather complex process, which has so far defied 
most attempts for a complete analytic treatment even for 
simple theories. It is a far-from-equilibrium phenomenon 
involving densities inversely proportional to the coupling. 
The nonperturbatively large occupation numbers cannot 
be described by standard kinetic descriptions. So far, 
classical statistical field theory simulations on the lattice 
have been the only quantitative approach available 
These are valid for not too late times, before the approach 
to quantum thermal equilibrium sets in. Up to now, stud- 
ies in quantum field theory have been mainly limited to 
linear or mean-field type approximations (leading-order 
in large- TV, Hartree) jfj, which present a valid description 
for sufficiently early times. However, they are known to 
fail to describe thcrmalization and miss important rescat- 
tering effects 0, • Going beyond mean- field has long 
been a major difficulty in practice: Similar to pertur- 
bation theory, standard approximations such as based 
on 1/iV-expansions of the one-particle irreducible (1PI) 
effective action can be secular in time and do not pro- 
vide a valid description. In contrast, it has recently been 
demonstrated 0, Q for 1 + 1 dimensional theories that 
far-from-equilibrium dynamics and subsequent thermal- 
ization can be described using a 1/iV-expansion of the 
two-particle irreducible (2PI) effective action [{j. Hoi Hl| . 
Below, we show that this provides a systematic and prac- 
ticable nonperturbative approach for quantum field theo- 
ries in 3 + 1 dimensions and with a nonzero macroscopic 
field (j>, relevant for realistic particle physics applications. 

In this work, we present the first quantum field theo- 
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retical description of the phenomenon of parametric res- 
onance taking into account rescattering: For an O(N)- 
symmetric scalar field theory we employ the 2PI 1 /TV- 
expansion to next-to-leading order (NLO), which in- 
cludes off-shell and memory effects 0, El- We point 
out that the approach solves the problem of an explicit 
description of the dynamics of correlation functions at 
nonperturbatively large densities. We present a complete 
numerical solution of the corresponding equations of mo- 
tion. Moreover, we identify the relevant contributions to 
the dynamics at various times and provide an approxi- 
mate analytic description of the nonlinear dynamics for 
the entire amplification range. Apart from the resonant 
amplification in the linear regime, we identify two charac- 
teristic time scales, which signal strongly enhanced parti- 
cle production in a broad momentum range due to nonlin- 
ear, source effects. This collective amplification is crucial 
for the rapid approach to a subsequent, quasistationary 
regime, where direct scattering drive a very slow evolu- 
tion towards thermal equilibrium. We emphasize that 
these processes cannot be seen in mean-field approxima- 
tions. Similar phenomena have been observed in classical 
field theories |5j , and we present an analytic criterion for 
the validity of classical statistical approximations. These 
effects are important for a reliable description of the sys- 
tem at the end of the resonance stage for finite N < 1/A. 
For realistic inflationary models with typically A <C 1 this 
is, in particular, crucial to determine whether there are 
any radiatively restored symmetries. 

We consider a relativistic real scalar field (p a (a = 
l,...,N) with action S[ip] = - j x {\ip a (U x + m 2 )ip a + 

4^7 (^Pafa) 2 }, where summation over repeated indices is 
implied. We use the notation J = J c dx° J dx with 
C denoting a closed time path along the real axis. All 
correlation functions of the quantum theory can be ob- 
tained from the 2PI generating functional for Green's 
functions T[4>, G], parametrized by the field expectation 
value <fi a (x) — ((p a (x)) and the connected propagator 
Gab{x,y) = (T c <Pa(x)<Pb(y)) ~ (j) a {x)4>b{y) E3 : 

m G] = S[<l>] + ^TrlnG- 1 + ^Tr G^) G + T 2 [0, G] , 
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FIG. 1: The dots indicate that each diagram is obtained 
from the previous one by adding another "rung" with two 
full propagator lines at each vertex. The crosses denote field 
insertions. 

where tG^foy;^) = ^XwUr The term T ^ G ^ 
contains all contributions beyond one-loop order and can 
be represented as a sum over closed 2PI graphs ^(J. To 
NLO in the 2PI 1/iV-expansion, r2[</>, G] contains the 
diagrams with topology shown in Fig. ^ ii H]- The 
infinite series can be summed analytically: 

r 2 [0,G] = -^ j G aa (x,x)G bb (x,x) + l -TrlnB(G) 



i\ 2 



(6^) 2 J xy 



B 1 (x, z; G)G 2 (z, y)cj) a (x)G ab (x, y)<j>b(y) 



with B(x, y; G) = S(x -y) + §?G 2 (x, y) and G 2 (x, y) 
G a b(x,y)G a b(x,y). The equations of motion |14| are 



5<p a (x) 







ST[(/),G} 
8G ab {x,y) 



0. 



(1) 



As we will show below, for the phenomenon of paramet- 
ric resonance, each diagram of the infinite series in Fig.^ 
eventually contributes at the same order in the coupling 
A. In this sense, the 2PI 1/iV-expansion to NLO rep- 
resents a minimal approach for a controlled description 
including rescattering. This justifies the rather involved 
complexity of the approximation. 

Overview. Parametric amplification of quantum fluc- 
tuations can be best studied in a weakly coupled 
system that is initially in a pure quantum state, 
characterized by a large "classical" field amplitude 
<i>a (t) — cr(t)Mo yj 6N/ A 5 a i, and small quantum fluctua- 
tions, corresponding to vanishing particle numbers at ini- 
tial time. Here Mo sets our unit of mass and the rescaled 
field a(t — 0) = (To is of order unity. We first present 
results from a full numerical |15| integration of the time 
evolution at NLO. In Fig. [5] we show the classical-field 
and the fluctuation contributions to the (conserved) to- 
tal energy as functions of time. The former dominates 
at early times and one observes that more and more en- 
ergy is converted into fluctuations as the system evolves. 
A characteristic time - denoted as t n onpert in Fig- HI 
is reached when the classical and the fluctuation part of 
the energy become of the same size. Before this time, 
the coherent oscillations of the field cj) lead to a reso- 
nant enhancement of fluctuations in a narrow range of 
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FIG. 2: Total energy (solid line) and classical-field energy 
(dotted line) as a function of time for A = 10~ 6 . The dashed 
line represents the fluctuation part, showing a transition from 
a classical-field to a fluctuation dominated regime. 

momenta around a specific value |p| — Pq: this is para- 
metric resonance 0, |g ■ Nonlinear interactions between 
field modes then cause this amplification to spread to 
a broad range of momenta. We find that the resulting 
rate of amplification exceeds the characteristic rate 70 
for the resonant growth. This is illustrated in Figs. |3| 
and 01 where the effective particle numbers Q are dis- 
played for various momenta as a function of time, both 
for the transverse (G±) and the longitudinal (Gy) sector, 
with G ab — diag{G||, G±, . . . , G±}. One observes that, 
in contrast to the rapid changes for particle numbers be- 
fore inonpert) a comparably slow quasistationary evolu- 
tion takes place at later times. The fluctuation domi- 
nated regime for t > i non pert is characterized by strong 
nonlinearities. For instance, from Fig. one infers for 
t ~ inonpert that the classical field decay "overshoots" 
and is temporarily reversed by feed-back from the modes. 
This can be directly seen in the evolution for the rescaled 
field shown in Fig. [5j The particle numbers of Figs. |3| 
and 01 exhibit correspondingly a reverse behavior. Very 
similar phenomena have been observed in related classi- 
cal theories [f| , though a direct comparison requires sim- 
ulations for the same model as studied here |l2(. There 
the full nonlinearities, i.e. including all orders in 1/N, 
are taken into account while leaving out quantum cor- 
rections. This indicates the capability of the 2PI 1/N 
expansion at NLO to capture the dominant nonlinear 
dynamics. 

The characteristic properties described above can be 
understood analytically from the evolution equations for 
the Fourier modes of the one- and two-point functions. 
Separating real and imaginary parts with Gu t ±(t, t'; p) = 

P) - h /°II,-l(*> P) si S n c(*-*')> the real P\\-X de- 
note the spectral and Fi\± the statistical two-point func- 
tions Ull. We define: 



M 2 (t) 



A 



z + —[3mt) + (N-l)T^t)], (2) 

where T\\ t± (t) = / A {g^ri J j|,x(t,*;p) with A>p . Ini- 
tially, F||(0,0;p) = l/2w||(p), 9 t F||(t,0;p)| t=0 = 0, 
dtdt>F\\(t,t';p)\t=t'=Q — w ||(p)/2, and similarly for F±. 
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FIG. 3: Effective particle number density for the transverse 
modes as a function of time for various momenta p < 5po. 
At early times, modes with p ~ po are exponentially am- 
plified with a rate 270. Due to nonlinearities, one observes 
subsequently an enhanced growth with rate 670 for a broad 
momentum range. 



The frequencies arewii(p) = [p 2 



-M 2 (1- 



Ml 



Af 2 (0), and similarly for cjj_(p) with 3er 2 — > a 2 . The ini- 
tial conditions for the spectral functions are fixed by the 
equal-time commutation relations: p\\±(t, t'\ p)\t=r = 

fmddtp\\^t,t'\v)\t=f = 1 Hi- 

(7) Early-time linear (Lame) regime: parametric reso- 
nance. At early times the cr-field evolution equation re- 
ceives the dominant (0(A )) contributions from the clas- 
sical action S. The classical field dynamics is character- 
ized by rapid oscillations with constant amplitude and 
period 2tt/ujo- The evolution equations for the two-point 
functions at 0(A°) correspond to free-field equations with 
the addition of a time dependent mass term 3cr 2 (i) for the 
longitudinal and a 2 (t) for the transverse modes. The dy- 
namics in this linear regime has been extensively studied 
in the literature and is known to be described by the so- 
lution of a Lame equation Q Q: Parametric resonance 
manifests itself by an exponential growth of the statisti- 
cal two-point functions describing particle production in 
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a narrow momentum range with p 2 < . Averaging over 
the short-time scale ~ Uq 1 one finds for the transverse 
modes for t, t' » ^q 1 : 



F ± {t,t';p ) 



= 70(*+O 



(3) 



for the maximally enhanced mode with po — ^f- O nc 
finds a much smaller growth for the longitudinal modes 
(cf. also Fig. 0J. The analytic solution to O(X ) agrees 
accurately with the NLO numerical results at early times. 

(II) Source-induced amplification regime: enhanced 
particle production for longitudinal modes. Because of 
the exponential growth of the transverse fluctuations for 
7o< 3> 1, the 0(A°) approximation eventually breaks 
down at some time. Using the two-loop graphs of Fig. ^ 
at 0(A) the evolution equation for Fji can be approxi- 
mated by 



(a 2 + p 2 + 3a 2 (t) + M 2 {t)) F,| (t, f ; p) ~ (4) 
2aa 2 {t) T ± (t) F]| (t, t'; p) + aa{t)cr(t r ) U.^(t, t'\ p) , 
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FIG. 4: Same as in Fig. [3] for the longitudinal modes. Non- 
linear source effects trigger an exponential growth with rate 
470 for p < 2po- The thick line corresponds to a mode in 
the parametric resonance band, and the long-dashed line for 
a similar one outside the band. The resonant amplification is 
quickly dominated by source-induced particle production. 



= J dqF ± (t,t';p 
with constant c 



q)Fj_(i,f';q), and 
1. Here, we have 



with n£(M';p) 

_ A N-l c 
a (2tt)3 6N Z% 

used that the momentum integrals are dominated by the 
amplified transverse modes. In addition, we exploit the 
fact that, for t" ~ t < t„ 



knonpert j 
?2 1, ,11. „ \ „ 2 



Fl(t ) t";p )>pi(t,t";p )/4. 



(5) 



Its validity can be seen from @ and the corresponding 
behavior of p ± (t,t";p ) ~ e 7o(t-t") f or t " < t. We em- 
phasize that strictly neglecting in the evolution equations 
/9 2 -terms as compared to F 2 -terms for all modes and all 
times corresponds to the classical statistical field theory 
limit This provides an explicit analytic criterion for 
the applicability of the classical methods employed in [|| . 
For sufficiently early times, all momentum integrals are 
indeed dominated by the enhanced modes p ~ po for 
which (0 is valid. A detailed analysis reveals that the 
memory integrals appearing in the equations of motion 
are dominated by the latest times t" ~ t: One can ex- 
plicitely verify that effective locality, described by the 
replacement J* dt" — > Jt_ c / Uo d*" Ec h ©)' is va ^ 
for the time-averaged behavior over the short-time scale 
~ Uq . We conclude that classical statistical approxima- 
tions can provide a good description in this regime. 

The first term on the RHS of 10} is a NLO contribution 
to the mass, whereas the second term represents a source- 
term, resulting from annihilations of amplified transverse 
modes as well as stimulated emission processes. To eval- 
uate the momentum integrals, we use a saddle point ap- 
proximation around p ~ po, valid for t, t' 3> , with 
Fi(t, t>, p) ~ F ± (t, f , po) exp[-| 7o '|(i + t'){p - p ) 2 /2]: 



T±(t) 



n£(f,f ; o) 



p 2 F ± (t,t;p ) 

2(7^7^)1/2 ' 

p 2 F 2 (t,t';p ) 

^\i' a '\(t + t')Y/ 2 



(6) 
(7) 



Here we wrote the source term for its maximum at p = 0. 
More precisely, it affects all modes with p < 2po. An im- 
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FIG. 5: The rescaled field a as a function of time for A = 10 
(up) and A = 10 (below). 

portant characteristic time is reached when the mass cor- 
rections become comparable to the classical mass term. 
This can be seen to happen at t no npert — ( m A _1 )/(27o), 
when Tj_ ~ 0(X~ 1 ) (see also Note that at this time, 
both the LO and the NLO mass terms are of the same 
order in A, but with opposite sign. We point out that the 
source term in Q becomes important at the earlier time 

^source — ^nonpcrt/2 5 (8) 

when F ± (t,t';p ) ~ O^X- 1 ' 2 ). For t > t sourcc , the 
exponentially growing source drives the dynamics of lon- 
gitudinal modes and one finds a strong amplification for 
p < 2p : Fu(t,t';p) ~ exp[27 (i + t')]. This agrees pre- 
cisely with the numerical results shown in Fig. 0] 

(III) Collective amplification regime: explosive particle 
production in a broad momentum range. A similar anal- 
ysis can be made for the transverse fluctuations. The ap- 
proximate evolution equation for Fx has a similar struc- 
ture as (@J. Beyond the O(X ) (Lame) description, it 
receives contributions from the feed-back of the longitu- 
dinal modes at 0(A) and from the amplified transverses 
mode at 0(X 2 ). The corresponding mass corrections re- 
main small until inonpert, whereas the source term is para- 
metrically of the form ~ X 2 F^_/N. This leads to the 
characteristic time 

^collect ^ 2t nonpcrt /3 + (lniV)/(67o) (9) 

at which F±(t,t';p ) ~ ©(iV 1 ^- 2 / 3 ). Correspondingly, 
for tcoiiect t < inonpert one finds a large particle produc- 
tion rate ~ 670 for a wide range of momenta, in agreement 
with the full NLO results in Fig. In this time range 
the longitudinal modes exhibit an enhanced amplifica- 
tion as well (cf . Fig. 0} . The abundant particle produc- 
tion is accompanied by an exponential decrease of the 
classical-field energy. We emphasize that the collective 
amplification regime is absent in the LO large- N approx- 
imation. Consequently, even for the transverse sector the 
latter does not give an accurate description at interme- 
diate times if i C ollect < inonpert, that is for N < A -1 . 

(IV) Nonperturbative regime: quasistationary evolu- 
tion. At t ~ inonport one finds F±(t,t';p ) ~ O(N A" 1 ). 



As a consequence, there are leading contributions (~ A ) 
to the dynamics coming from all loop orders (cf. Fig.^l. 
In particular, the evolution equations are no longer "lo- 
cal" in the sense described under (II) and memory ef- 
fects become important. In contrast to the rapid reso- 
nant dynamics before tnonpert, a comparably slow, qua- 
sistationary evolution driven by direct scattering sets in. 
We emphasize that the collective amplification regime 
triggered a rapid approach to monotonously decreasing 
particle number distributions as functions of momentum, 
which become quasistationary afterwards (cf. Figs. El II}- 
The approach to true thermal equilibrium is exceedingly 
slow for the employed range of couplings A = 1CP 6 — 10, 
for which a non-negligible parametric resonance regime 
can be observed. For phenomenological applications it is 
therefore crucial that in the quasistationary, pre-thermal 
regime in the spirit of Refs. [2U121] a number of aspects do 
not differ much from the late-time thermal regime [l3| . 
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